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Abstract 

We study the dynamics of an elastic body whose shape and position 
evolve due to the gravitational forces exerted by a pointlike planet. The 
main result is that, if all the deformations of the satellite dissipate some 
energy, then under a suitable nondegeneracy condition there are only three 
possible outcomes for the dynamics: (i) the orbit of the satellite is un- 
(~| , bounded, (ii) the satellite falls on the planet, (iii) the satellite is captured 

^ f in synchronous resonance i.e. its orbit is asymptotic to a motion in which 

the barycenter moves on a circular orbit, and the satellite moves rigidly, 
always showing the same face to the planet. The result is obtained by 
making use of LaSalle's invariance principle and by a careful kinematic 
analysis showing that energy stops dissipating only on synchronous orbits. 
We also use in quite an extensive way the fact that conservative elastody- 
namics is a Hamiltonian system invariant under the action of the rotation 
group. 

1 Introduction 

In this paper, we study the dynamics of an elastic satellite interacting with a 
pointlike planet. Precisely, we study the dynamics of an elastic body, moving in 
the gravitational field generated by a pointlike mass. We consider the equations 
of motion of continuum mechanics, with body forces due to the gravitational 
fields and internal traction arising from the body deformation, without intro- 
ducing any further approximation. We prove that, if the internal structure of 
the satellite is such that any deformation dissipates some energy and if a suit- 
able nondegeneracy condition is satisfied, then the dynamics of the system has 
only three possible final behaviors: 

(i) the orbit of the satellite is unbounded; 

(ii) the satellite falls on the planet; 

(iii) the satellite is captured in synchronous resonance. 

By item (iii) we mean that the shape of the body reaches a final configuration, 
that its center of mass moves on a circular orbit and that it always shows the 
same face to the planet, i.e. the planet is at rest in a frame comoving with the 
satellite. 
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Concerning the inner structure of the body we make as few assumptions as 
possible. Precisely we assume that the stress tensor is the sum of two terms, 
the first one being Hamiltonian, i.e. it is the gradient of a "stored energy 
functional", and a second one being nonconservative. On the second term we 
only assume that, as a consequence of its presence, there is dissipation of energy 
at any time at which the time derivative of the Cauchy Green stress tensor does 
not vanish. 

The idea of the proof is to use LaSalle's principle (see I15j), which is a gener- 
alization of Lyapunov theorem. LaSalle's principle ensures that any precompact 
orbit approaches an invariant set which is contained in the manifold where the 
Lie derivative of the energy vanishes. The core of the paper consists in char- 
acterizing such an invariant set. Since in such an invariant set the dynamics is 
conservative, it turns out that a convenient framework for our study is that of 
Hamiltonian systems with symmetry as developed for example in [18j or, in a 
form directly useful for our problem, in [20j . 

So we start by writing down the Lagrangian and the Hamiltonian of the 
conservative part of the system and then we add to the equations of motions 
the nonconservative forces. 

Then we start analyzing the nondissipating manifold MV. We first prove 
that consists of rigid motions, and then we show that the motions laying 
on are actually circular orbits. Finally we show that they are relative equi- 
libria of the reduced Hamiltonian system obtained by exploiting the rotational 
invariance of the original Hamiltonian. At this point the application of LaSalle's 
principle would allow to conclude that the orbit is asymptotic to a manifold ob- 
tained by taking the union of all the synchronous orbits. Li order to prove that 
the system is actually asymptotic to a single synchronous orbit we exploit the 
conservation of angular momentum and we assume a nondegeneracy property 
stating that the relative equilibria are isolated. This property is discussed in 
detail in Section S] and we show that it is typically fulfilled. 

The present result still has some quite strong limitations. The main one is 
that we do not discuss existence and uniqueness for the Cauchy problem of the 
equations we study. As is well known this is one of the main open problems of 
elasticity, so we do not enter its discussion. Here we limit ourselves to assuming 
that the system we study is well posed. Our assumptions are true for example for 
any Galerkin cutoff of the system or if one can show that the system behaves as 
a parabolic system (as is expected in the dissipative case under consideration). 

The second limitation of our result rests in the fact that we assume that the 
system is described by differential equations. This means that we do not consider 
the case where the system is described by an integrodifferential equation with 
delay, a case which can occur in elasticity. The case we have in mind is the 
one in which the dissipation is of the kind of that appearing in Navier Stokes 
equation. We expect that our theory can be extended to the case with delay, 
but for sure the methods should be adapted. 

The study of the gravitational interaction between a deformable body and 
a pointlike mass traces its origin back to the pioneering work by Darwin [7,8]. 
His work shows that, in some approximation, the effect of the internal dynamics 
of the satellite is just that of producing an effective dissipating effective force 
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on the orbital and spin degrees of freedom of the satelhte. Darwin's work was 
subsequently generalized by Kaula |13| and many other authors (for instance, 
[21[T21[I71[Tnj ) . Critical reviews of the work by Darwin, Kaula and followers can 
be found in However, the Darwin-Kaula procedure is heuristic and, from 

a mathematical point of view, its range of validity is far from being clear. For 
this reason, in the present paper (as in [Bl), the point of view is that of starting 
from first principles in order to obtain a rigorous mathematical proof of the 
phenomena under consideration. 

We remark that the result of the present paper rules out the possibility 
that periodic orbits different from synchronous resonance exist. This is quite 
surprising, since some celestial bodies are known to be in a spin-orbit resonances 
different from the synchronous one (for example Mercury). Our interpretation 
is that orbits like the 3:2 resonance of Mercury are very likely to be either 
metastable or orbits eventually impacting the planet. We plan to investigating 
further in this direction. 

In Section [2] we state the main result of the present paper: to this end, we 
recall the Lagrangian formalism for elastodynamics, we write down the related 
Cauchy problem and we formulate the nondegeneracy assumption. Section [3] is 
devoted to the proof of the main result: we recall the statement of LaSalle's 
invariance principle, prove that the only solutions which dissipate no energy 
are synchronous orbits and apply La Salle principle to our system. Finally, in 
Section |4] we discuss the nondegeneracy assumption and we prove that typically 
it is fulfilled. 

2 Statement of the main result 
2.1 The setting 

We study the dynamical system consisting of: 

(i) a pointlike mass M (which we will sometimes call "planet"), which is at 
rest and which is chosen as the origin of a system of coordinates; 

(ii) an elastic body, free to move in space; we will call this extended body 
"satelhte". 

To deal with elastodynamics we use the framework of [18] and [20J from 
which we take some notations and formalism, that we now recall. 

We denote by S C the reference configuration of an elastic body and 
assume that B is open and bounded with a smooth boundary dB. We define the 
configuration space Q to be a Banach space of map^ : B ^ 'R^. Typically in 
elastodynamics one assumes Q C H'^ {B) with s large enough; we will come back 
later to this point, for the moment we simply assume that C, admits as many 
derivatives as needed. 

In the conservative case, classical three dimensional elasticity is a Lagrangian 
system, the Lagrangian £ : TQ — > 3R is the difference of kinetic and potential 

^Actually we should restrict ourselves to the manifold of the maps s.t. det DC, > 0, however 
we will consider this as a condition on the domain of definition of the system. 
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energy. In our case there are also some dissipative forces that will be added to 
the Lagrange equations. As usual TQ ~ Q © Q is the tangent bundle to Q. 

We start by writing down the conservative part of the system. The La- 
grangian C of the system is defined by 



where 



C = K -Ug-Usg-Ue , (2.1) 

K{0--^1 I Po(x) C(x)'d^x (2.2) 
UgiC) / po(x)K,(C(x))d3x , (2.3) 

JB 

UsAO / Po(x)V^4(C(x))d3x , (2.4) 

JB 

C/e(i^C) / VK(x,i?C(x))d'x , (2.5) 

JB 

the functions Vg, V^'g are defined by 

Vgix) ■■= , (2.6) 

and W is the stored energy function; k is the universal gravitational constant, 
and po € C°°{B) the density of the body in the reference configuratior^. The 
stored energy function is assumed to depend on ^ only through the deformation 
gradient F := Z?C = {dC /dx"^}. We assume that W is frame independent in 
the sense that 

VF(x, F) = Wix, RF) for aU R e 50(3) . (2.8) 
As shown in [20J this implies that the Kirclioff stress tensor, namely 

dC dW 

" d^didOjd^ ^''''^ " uiiderstood) , (2.9) 

is symmetric. 

Assuming the stress free boundary condition, namely 



dW 



= , (2.10) 

dB 



where n = (ni, 712,713) is the external normal to dB, one deduces the standard 
Lagrange equations: 

PoC = ~VcC ^ ~po^(C) - Po^(C) + Q^Q^Q^^ - (2-11) 

where V^i2 = (Vri£, V^2£, V^3£) is as usual the gradient with respect to the 
scalar producl|f| and coincides with the expression at r.h.s. of (I2.1ip . 

^Of course we assume that po(x) ^ Vx g B. 

3i.e. it is defined by dC{C)h = (V£{C); /i>i2 for all h e Q. 
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Remark 2.1. It is easy to check that, if a function C/(C) is rotation invariant, 
I.e. U{RC) = U{C), Vi? e SO{3), then 



[VU]{RO = i?Vf/(C) . (2.12) 
All the terms of the Lagrangian have this property. 

Since the Lagrangian is independent of time, the energy 

n^K + Ug + U,g + Ue . (2.13) 

is formally conserved for the system (j2.1ip with the boundary conditions (|2.10D . 
Furthermore, since the Lagrangian is invariant under the group action 

Q X 5(9(3) 3 (C, i?) H> i?C, (2.14) 
by Nother's theorem the quantity 

L := / C(x) X po(x)C(x)d3x , (2.15) 

JB 

is conserved for the system. Of course L coincides with the total angular mo- 
mentum. 

In order to get the equations governing the non conservative dynamics one 
has simply to add the nonconservative forcesjfl, G = G(C,Ct) i-S- to substitute 
equation (|2.1ip with the equation 

PoC - -Vc/: - G . (2.16) 

In order to write down the precise assumptions on G (which will be given 
in the next subsection) we have also to introduce the (right) Cauchy Green 
deformation tensor C :— {DC,)^ DC^ or, componentwise 

^ dx'' dx'' 

a 

As it is well known C is symmetric, positive definite and allows to write DC, in 
the polar decomposition form DC,{x.) — i?(x)-y/C(x), where i?(x) G 50(3) is a 
rotation matrix. 



2.2 The Cauchy problem 

In the following we will always denote by y = (C, Q G TQ ~ Q ® Q a point in 
the space of initial data for our system (which we keep distinct from the phase 
space, in which the velocities will be substituted by the momenta). 

The problem of existence of solutions for the system (j2.11l) . (|2.16p with 
the boundary conditions (I2.10p is a major problem of elastodynamics (see e.g. 
[5]). Here, we do not want to enter such problem, so we will limit ourselves 
to assuming the needed well-posedness properties. Furthermore we will only 
consider classical solutions, so we give the following definition: 

*by this notation we mean that G is a function of the functions (^t, not of their value 
<^(x), CtC'^): so it can also depend on an arbitrary number of derivatives of such functions. 
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Definition 2.2. Givenyo G TQ, apositiveT and a function y G C^((0,T);TQ) 
we say that it is a solution of the system (j2.16p . (|2.10p with initial datum yo, 
if it fulfills the equations and the boundary conditions for all t € (0, T) and one 
has 

lim y{t) = yo . 

Then we also need the following definitions: 

Definition 2.3. A solution is said to be impacting (in the future) i/inft>o dist(^(;B, t), 0) = 
0. 

Definition 2.4. A configuration is said to be non singular if det[D(^] > 0. 

Definition 2.5. A solution y{t), t e (OjT) is said to be regular if it is non 
impacting and the corresponding configuration is non singular for all times. An 
initial datum yo = (CoiCo) regular if dist{(^Q{B),0) > and det[D(^o\ > 0. 

Definition 2.6. A regular solution y{t), t e (0,oo) is said to be precompact if, 
for any increasing sequence {t,i} C (0, oo) there exists a subsequence {^rifc} s.t. 
the limit limfe_j.+oo 2/(^rifc) exists and the limit is regular. 

Assumption 1. We assume that 

(i) for all regular initial data yo £ TQ, the Cauchy problem for the system 
(|2.16D with the boundary conditions (|2.10p is locally well-posed; 

(ii) let y{t) be a non-impacting solution. Then its time of existence is infinite 
and it is forever non singular. 

(Hi) Any regular solution fulfilling sup^^g |C(a;, t)| < oo is precompact; 

(iv) the angular momentum L is conserved along the solutions; the Lie deriva- 
tive of the energy (|2.13p is nonpositive and vanishes if and only if C = 0, 
where C is the Cauchy Green tensor (cf. eg. (j2.17p j. 

Remark 2.7. If one thinks that the nonconservative equations we are studying 
behave like parabolic equations, then typically the situation of Assumption [7] 
is fulfilled. Indeed, for positive times the solution of a parabolic equation is 
contained in all Sobolev spaces. Existence of global solutions has been proved 
for some equations of viscoelasticity (see JJ^ and flW). Moreover, the property 
of having a compact attractor is known for many classes of evolution equations 
involving dissipation (see J23^). 

Remark 2.8. A situation in which Assumption]^ is fulfilled is that in which 
the space Q is finite dimensional. A typical situation we have in mind is that in 
which the space Q is composed by maps obtained by Galerkin cutoff from some 
of the maps belonging to the original infinite dimensional configuration space. 
For example one could decide to keep only a finite ( arbitrarily large ) number of 
spherical harmonics of the maps describing the configuration. 
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2.3 The nondegeneracy assumption 

In the following (see subsect. 13. 2p we will prove that the nondissipating orbits 
are relative equilibria of the Hamiltonian system obtained by Legendre trans- 
forming the Lagrangian (j2.ip . We are now going to recall the notion of relative 
equilibrium and to state the nondegeneracy condition we need. 
Define the momentum 

TTi-poC, (2.18) 

then the Hamiltonian of the system coincides with the function "H given by 
(|2.13p . where however 

K - I Ij^'^^x ; (2.19) 

and the Lagrange equations (j2.11l) are equivalent to the Hamilton equations of 

The Hamiltonian is invariant under the action of the symmetry group SO{i) 
defined by 

Rz = R{tt, C) := (Rtt, RC) , (2.20) 

and the total angular momentum L (written in terms of positions and momenta) 
is the corresponding conserved quantity. Then one can use Marsden-Weinstein 
reduction procedure, that can be summarized as follows. 

(1) Fix a value Lg of L and consider the manifold 

A^Lo :={z = (7r,C) : L(z) = Lq} ; 



(2) Consider the subgroup C/Lq ^ 50(3) leaving invariant A^Loi namely the 
group of the rotations around the axis Lq. Consider the quotient mani- 
fold A^Lo/^/Lo- Such a manifold has a natural symplectic structure and 
furthermore the Hamiltonian T-L, as well as its Hamilton equations, pass 
to the quotient and define a Hamiltonian system on A4-l„/G'Lq- 

We denote by the Hamiltonian of such a reduced system. 

Definition 2.9. The critical points ofN-Lg are called relative equilibria of the 
Hamiltonian system %, at angular momentum Lq. 

Definition 2.10. A relative equilibrium is said to be topologically nondegenerate 
if it is not an accumulation point of relative equilibria with the same angular 
momentum. 

By abuse of notation a representative of the equivalence class of a relative 
equilibrium is also called a relative equilibrium of %. 

Remark 2.11. It is well known (see e.g. \Tj) that z^ is a relative equilibrium 
if and only if the Hamiltonian vector field of H. at z,, is tangent to the orbit of 
SO(3) through z^. 

Remark 2.12. If z,, is a relative equilibrium then the corresponding orbit z{t) 
( under the flow of the Hamiltonian system H) is formed by relative equilibria. 
In the nondegenerate case there are no other relative equilibria with the same 
angular momentum in a neighborhood of the orbit Utz{t). 
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Of course a relative equilibrium Zg corresponding to a value Lq of L is topo- 
logically nondegenerate if and only if the same is true for the relative equilibrium 
Rze, where R G SO{3) is arbitrary. 

Definition 2.13. A value £ Cz ^ of the modulus of the angular momentum is 
said to be nondegenerate if all the relative equilibria with angular momentum L 
satisfying |L| =i are topologically nondegenerate relative equilibria ofH. 

Remark 2.14. In Section^we will comment on this condition and show that 
it is in general fulfilled. 

Theorem 2.15. Let y(t) be a solution of eg. (j2.16p with the boundary condi- 
tion (|2.10p , and let I be the corresponding value of the modulus of the angular 
momentum. Assume that i is nondegenerate, then, one of the following three 
(future) scenarios is possible: 

(i) the trajectory of B is unbounded; 

(ii) the solution impacts the planet; 

(Hi) the solution is asymptotic to a synchronous non- dissipating orbit, which is 
a relative equilibrium with angular momentum t. 

3 Proof of theorem 12.151 
3.1 LaSalle's invariance principle 

In order to study the dynamics of the system, we make use of the LaSalle's 
principle which is a refinement of the classical Lyapunov's theorem. We now 
recall its statement and proof. 

Let 3^ be a Banach space and let Z// C 3^ be open. Consider a system of 
differential equations 

y = /(y) y^U , (3.1) 
We denote by the flow of which we assume to be locally well defined. 

Definition 3.1. Let 7 be the orbit of (j3.ip with initial condition i/q. A point rj 
is said to be an uj-limit point of j if there exists a sequence of times i„ — > +00 
such that 

lim ^'-{yo)=v. (3.2) 

n— >-+oo 

Definition 3.2. The uj-limit set of an orbit 7 is defined as the union of all the 
uj-limit points ofj. 

Definition 3.3. A solution y[t) d U, t (z (0, 00) is said to be precompact if, for 
any increasing sequence C (0, 00) there exists a subsequence {tn^} s.t. the 
limit limfc^-|_oo yitn^) exists and belongs to lA. 

Remark 3.4. It is well known that the uj-limit of a precompact orbit is a con- 
nected set (see e.g. IBI). 

The classical version of LaSalle's principle may be stated as follows: 
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Theorem 3.5. Suppose that Ti : U ^ ^ is a real-valued smooth function, 
such that £ fHiij) < 0, Vy G U, where £ f is the Lie derivative. Let X he the 
largest invariant set contained in MT> := \y G U\£ fT-L{y) = 0}, then the uj-limit 
of every precompact orbit is a non-empty subset of I. 

Proof. Let 7 := {ip*{yo)\t > 0} be a precompact orbit, and let T C U be the 
w-limit of 7. We prove now that T is invariant. Indeed, let G F, then there 
exists a sequence t„ — ^ +00 such that ip^"{yQ) rj. But we have 

ip^T^) = ip'i lim if^-iyo)) = lim <^*+*"(yo) e T . 

We prove now that the w-hmit is contained in AfV. Let rjo eT. Then there 
exists a sequence tn +00 such that (/3*"(yo) Vo- Now, let 

c:=Hir,o)= 1™ ^[(^'"(yo)] • 

n— >+cxD 

Since 'H[(p^{yQ)] is a time-nonincreasing function, lini„_>+oo H[v3*" (yo)] = c is 
independent of the subsequence t„ , and thus Ji = c on the whole F. By the 
invariance of F it follows that = on F, and therefore 7 C I. □ 

3.2 Non-dissipating orbits 

Consider the non dissipating manifold defined by 

MV:^{yeQ(BQ : n{y) = o} , (3.3) 

where, for short we denoted by lily) the Lie derivative of H along the vector 
field corresponding to the equations (|2.16p . In this section we prove that the 
subset X C MT) invariant under the dynamics is formed by relative equilibria of 
the Hamiltonian system (|2.13p . 

Remark 3.6. On MT> the Lagrange equations (j2.1ip coincide with the non 
conservative equations (I2.16p . 

First we prove that the body is rigid along any orbit in I (we think that this 
should be well known, but we were not able to find a reference). 

Lemma 3.7. Let y G C^{{0,+oo),Q® Q) be a solution of ([^1^ s.t. y{t) = 
{C{t), C{t)) £ A/"!? yt G (0, +00), then, for all x,yeB, one has 

||C(a;)-C(y)l=o. (3.4) 

Proof. Fix two arbitrary points x, y G S and consider a path 7 C S connecting 
X to y. Let s be the arclength parameter, so that |^^^^| — 1 and in the reference 
configuration the path 7 has length 

£(7) ^ I ds. 

The length of the deformed path C(7) is expressed in terms of the Cauchy Green 
tensor C by (see e.g. [21J ) 

^(C(7)) = / 
J7 



ds 



C 



ds 



ds 



(3.5) 
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Therefore, since C = 0, we have ^[^(C(7))] = 0. Now, take two arbitrary times 
to,ti and let CojCi be the corresponding body configurations. We have that 
Ci ° (Co)"""^ : Co{S) Ci('S) is a length-preserving map between the Co{B) C K'^ 
and Ci{B) C 3?'^, both equipped with the restriction of the Euclidean metric on 
Moreover, o (Co)~^ is a diffeomorphism. 

Using the fact that the segments minimize the distance it is easy to conclude 
the proof of the lemma (some care is needed in order to take care of the fact 
that Ci (B) could fail to be convex) . □ 

Corollary 3.8. Let y = (C, C) ^■^ ^'^ statement of Lemma \3.7\ then there 
exist ^eQ, R€ C2((0, +oo), SO{3)) and Y e C^{{0, +oo), SR^) g.t. 



Cix^t) = Rmix) + Y{t)) . 



(3.6) 



The reference frame with origin Y(t) and coordinate axes R{t)ei is usually 
called comoving frame. In this frame all the points of the satellite are at rest 
along the orbit y{t). In particular — Y(i) is the position of the planet M in the 
comoving frame. 

Lemma 3.9. Let y G C'^{{0,+oo), Q® Q) be a solution of s.t. y{t) = 

(C(^)iC(O) G A/"!? Vi G (0, +cxd), then the quantity Y(t) in p.6p evolves in such 
a way that 



yij,k, 



:{Y{t))^G 



(3.7) 



is independent of time. 



Proof. Inserting the expression p.6p in the Lagrange equations (|2.1ip and ex- 
ploiting the rotational invariance of the r.h.s. (cf. Remark 12. ip one gets the 
following equation for ^, R and Y: 



Y +{Cj + ujuj)Y + 2wY +{Cj + lulo)^ 



1 d 



dW 



dx dx Po dx'^ d{d(^/dx^) 

where as usual a) := R^R. Denote for short 

L := {lj + idib) 



(DO 



(3. 



(3.9) 



and take the time derivative of 
depend on time one gets 



Taking into account that ^ does not 



dt 



LY + 2loY 



Li 



d 
It 



dx 



(C + Y) 



(3.10) 



Take now the derivative of such a quantity with respect to x", one gets the 
componentwise equation 



Ed 
It 

k 



dx'dx'' 



(3.11) 
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or using the invertibihty of the matrix 



a 



.(C(x)+Y(i)) 



(3.12) 



This equation shows in particular that the r.h.s. is independent of x. Due to 
the invertibihty of ^ and to the analyticity of Vg this means that the function 
of X 



d 
di 



(3.13) 



is actually independent of x- Thus taking the derivative with respect to and 
evaluating at x = one gets the thesis. □ 

Lemma 3.10. Let y e C^((0, +oo), Q® Q) be a solution of ((TTC)) s.t. y{t) = 
(C(OiC(i)) £ J^'D G (0, +oo), then the quantity Y(t) in (|3.6p is actually 
independent of time. 

Proof. We write down explicitly ^J}. We denote {Y'' {t) .Y"^ {t) {t)) = Y{t) 

and (Fl)2 + (r2)2 _^ (y3)2 ^ ^2^ Q^g j^^g 

d^Vg ^ ikMY\^{Y^f^ir^) ^^^^^^ 



d{x^Y r 
d'^Vg _ 3fcMr2(5(yi)2 ^2) 



dix'Ydx^ W - (3.15) 



9(x^)^i9x^ r'^ 



r 



Choose now the comoving frame in such a way that Y^ = Y^ ^ Q and 
at t = (which is possible up to redefinition of R and ^) . By (|3.15p and (j3.16p , 
for any t one must have 

Y'^(b{Y^ f ~r^) = 

F3(5((ri)2 -r^) =0 

which by continuity implies Y'^{t) = Y^{t) = and ^^(t) = r{t). Finally, 
substituting in (j3.14L one has 



{Y^{t)Y {Y^Y ' 

whose only solution is Y^{t) = ^o^. □ 

Remark now that, given a non dissipating solution, one can associate to it a 
shape of the body described by the function ^(x) + Y, and the shape evolves by 
a rigid motion about the fixed point M . Introduce the angular velocity which 
is defined as usual as the vector lo s.t. the two operators 



coincide. Then the velocity of the motion, in the comoving frame, is given by 
w X (C + Y). 

We have now that, for a non dissipating solution, w does not depend on time. 
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Lemma 3.11. The angular velocity io of a nondissipating solution y is inde- 
pendent of time. 

Proof. Consider again equation p.l2p . Since we now know that Y is indepen- 
dent of time it follows that the operator L (cf. eq. (j3.9p l fulfills L = 0. This 
means that, for any vector x one has 

^[ujxx + ujx{ujxx)]^0 . (3.17) 

To exploit such an equation take x — ^^nd project the square bracket on e^. 
Using standard vector identities this implies 

^ |w X e,;|^ = Vi . 

A trivial computation shows that this implies io = 0. □ 

Corollary 3.12. Let y{t) be a nondissipating solution as above, then it is the 
orbit of a relative equilibrium of the system (|2.13p . 

Proof. We have proved that along a non dissipating solution C,{t) = R{t)Co with 
a suitable configuration Co and a rotation matrix R{t) that we can choose in 
such a way that i?(0) = /. It follows C(t) = R{t)[uj x ^] = i?C(0). Passing to the 
phase space one gets that along such an orbit Tr{t) = po({t) — i?(t)7r(0). This 
shows that the solution is actually an orbit of the symmetry group, and this is 
a characterization of being a relative equilibrium. □ 

Thus we have that the manifold I is the union of the trajectories of all the 
possible relative equilibria of the system. 

3.3 End of the proof 

Applying La Salle principle to our system, with U defined to be the set of regular 
configurations, we get the following Lemma. 

Lemma 3.13. For any solution to (I2.16P with the boundary condition (j2.10p . 
one of the following three (future) scenarios occur: 

(i) the trajectory is unbounded; 

(ii) the solution impacts the planet; 

(Hi) the solution is asymptotic to the non- dissipating invariant manifold I. 

End of the proof of Theorem \2.15\ Let Lo be the initial value of the angular 
momentum, then a regular bounded solution is asymptotic to 

Xn{(C,C)eTQ : L(C,C)-Lo} , (3.18) 

but, by the nondegeneracy assumption the set (j3.18p is formed by orbits which 
are isolated in the invariant manifold |(C, C) £ TQ : L(C, C) = Lo|. Thus, by 
Remark 13.41 the w-limit set of an orbit is a single orbit in the set (I3.18p . i.e. a 
synchronous orbit. □ 
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4 On the nondegeneracy assumption 



In this section we are going to prove that if the restoring elastic forces described 
by the potential (|2.5p are strong enough, then the nondegeneracy of a relative 
equilibrium for the system (j2.ip is equivalent to the nondegeneracy of the rel- 
ative equilibrium for a rigid body having the shape given by the asymptotic 
configuration of the satellite. 

For simplicity, in this section we limit the discussion at the formal level, 
namely we forget all the difficulties related to the existence of unbounded oper- 
ators. All what follows is rigorous if Q is finite dimensional. It can also be made 
rigorous in the case of PDEs by detailing most of the assumptions, following 
the ideas of [5] and exploiting the ellipticity properties of the elasticity tensor 
(see [IB]), however this is outside the aims of the present paper. 

First of all, having fixed a configuration ^ we introduce the dynamical system 
describing the evolution of a rigid satellite with shape (. The configuration space 
SO{3) X SR'^ 9 (i?, x) and the dynamics is obtained from the Lagrangian obtained 
by restricting the Lagrangian (|2.ip to the set of motions of the form 

at)^R{t)[C + x{t)] ■ (4.1) 

Denote by £^ such a Lagrangian. One also has a corresponding Hamiltonian 
system, which is deduced from (12. ip in the standard way and whose Hamilto- 
nian coincides with the restriction of the Hamiltonian (j2.13p to the phase space 
of the rigid body. Denote by the Hamiltonian of the rigid body. Such a 
Hamiltonian is invariant under rotations so one can pass to the reduced system 
and to introduce again the relative equilibria and define the nondegenerate rel- 
ative equilibria for such a system and a nondegenerate value of the modulus of 
the angular momentum. 

In order to make a connection between the nondegeneracy of the relative 
equilibria for rigid motions and the nondegeneracy of the elastic motions we 
need to specify an assumption on the elastic potential energy Ue- Essentially 
we are going to assume that the elastic potential has a very steep, isolated (up 
to the symmetries) minimum at some shape. 

First remark that Ue is invariant under the action 

(50(3) X 5R3) X Q 9 {{R, x)X) ^ R[C + x] e Q , (4.2) 

and, given a point (, consider the group orbit C/^ := (5*0(3) x 5R'^)C C Q, so, if 
is a critical point of Ue, then all the orbit is critical for Ue- 

Assumption 2. One has Ue = ■^Ue, and Ue is a smooth function invariant 
under the group action (j4.2p with the further property that the set of its critical 
points is formed by finitely many orbits and each critical point is nonde- 

generate in the direction transversal to the group orbi^. 

Under this assumption we have the following 

Proposition 4.1. Fix a value Lq of the angular momentum, and assume e 
is small enough; let Q be a relative equilibrium of the Hamiltonian system H 

^Transversal nondegeneracy means that the restriction of Ue to any hyperplane transversal 
to the group orbit has a differential which is an isomorphism. 
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with angular momentum Lq. If (i?, 0) 6 5*0(3) x SR'^ is a nondegenerate relative 
equilibrium for the rigid system with Hamiltonian H^t then Q is a nondegenerate 
relative equilibrium for the elastic Hamiltonian system with Hamiltonian H . 

Proof. The proof is based on ideas from Lyapunov-Schmidt decomposition (see 
e.g. O). First of all it is useful to introduce suitable coordinates in Q in a 
neighborhood of (following the ideas of [6J). They are constructed as follows: 
let S C Q be a codimension 6 affine subspace transversal to the group orbit Q,^^ , 
then a suitable set of coordinates about Q is locally obtained by the map 

^xSO{-i)x^^3{^,R,x)^Ri^ + x) ■ (4.3) 

We now recall (and adapt to the present situation) some results of |20j . 
Using the method of Lagrange multipliers one immediately sees that the relative 
equilibria can be obtained by finding the critical points of — w • (L — Lq) under 
the condition 

L = Lo . (4.4) 

Here uj is the Lagrange multiplier. In [20J it was shown that this is equivalent 
to finding the critical points of the "augmented Hamiltonian'' Hi^„ defined by 

i^Lo := ^Lo + Uu, , (4-5) 

^ Jb Po 
C/l„ :=C/-i^Po|c^xC|'rf'x (4.7) 

again under the condition (|4.4I) (actually in the present case this is a straight- 
forward computation) . As pointed out in [20] the interest of this formulation is 
that the equations for the critical points of Hi,^ take the form 

TTe = paU} X Ce , (4.8) 
VC/(Ce) + POOJ X (W X Ce) = . (4.9) 

In particular the second equation is independent of tt. This allows to study 
separately (|4.9p . To this end we use the set of coordinates (|4.3p . However one 
has to pay attention to the fact that in general the section E is not orthogonal 
to the group orbit, so first we rewrite (|4.9p in the (original) dual form: 

dU{Ce)h+{poux{ux<:e);h)^2=0 , VheQ, (4.10) 

which in terms of the coordinates (|4.3p takes the form (at R = I) 

d^Ug{^e+Xe)h^ + {PQUJx{ujx{i^+Xe);h^)^^=Q , Vh^eU\ (4.11) 

d^UgiCe + Xe) + Usg{^e + Xe) + -^Ue{^e + Xe^^ (4.12) 
+ {POCJ X (CJ X (Ce + Xe)); h^) = , V/lJ S T^, E . 

which have to be solved together with 

TTe = PoUJ X i^e+Xe) , (4.13) 
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and the condition (|4.4p . In particular the system (|4.1ip . (|4.13p . (|4.4p is identical 

to the system for the reduced equilibrium of the rigid body with shape ^e, so 
by Assumption [2] it determines uniquely (up to a finite choice) Xe and lo (and 
TTe). We analyze now (|4.12p . Of course it is a perturbation oi d^Ue — 0, whose 
critical points are all nondegenerate (as functions of ^ they are nondegenerate in 
the standard sense), so, by the implicit function theorem a solution of (|4.12p 
must be close to a critical point of Ue and is also nondegenerate. This concludes 
the proof of the proposition. □ 

So we have reduced the problem of checking nondegeneracy to the problem of 
checking nondegeneracy of the relative equilibria for the motions of rigid bodies. 
This problem has been studied for example in |22j who obtained a complete 
characterization of the relative equilibria of a triaxial rigid body, provided the 
gravitational potential is approximated by its quadrupole expansion. In [22] 
the author obtained that, provided the distance Xe of the center of mass from 
the planet is large enough, there are exactly 24 families of stationary points of 
the reduced system. These stationary points are such that the principal axes of 
inertia are one pointing to M and a second one in the plane orthogonal to the 
plane of motion. The number 24 appears as the number of possible choices of 
the orientations of the body with prescribed principal axes of inertia. 

Remark that in particular it turns out that such critical points are nonde- 
generate. Furthermore we expect that, if Xe is large enough, then it should be 
possible to use the implicit function theorem to prove that the critical points 
are nondegenerate also for the system in which the gravitational potential is not 
subjected to any approximation. 

Remark 4.2. In 161 a problem is studied in which Assumption\^is violated due 
to the fact that the satellite is assumed to have a spherically invariant reference 
equilibrium configuration. In that case it was shown that actually the orbit is 
not asymptotic to a single synchronous orbit; nevertheless, one can still say that 
also in that case the orbit is asymptotic to a synchronous resonance, but the 
principal axes of inertia could rotate in the satellite. 
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